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A MODIFIED NEWTON-RAPHSON ANALYSIS OF FLIGHT MEASUREMENTS 
OF THE TRAILING VORTICES OF A HEAVY JET TRANSPORT 

By Lawrence W. Taylor, Jr., and Milton D. McLaughlin 
Langley Research Center 

SUMMARY 

A modified Newton-Raphson method has been used to determine the parameters in 
equations describing the vortex flow to obtain a best match with flight measurements of 
the flow behind the C-5A airplane. The flight measurements were made using a specially 
instrumented T-33 airplane which passed as closely as possible to the centers of the trail- 
ing vortices at several distances behind the C-5A airplane. The flight measurements were 
transformed to flow velocity relative to an inertial frame of reference. The assumed form 
of the flow consisted of the superposition of two counterrotating, finite core vortices. The 
positions of the vortex centers, their total circulation, the effective eddy viscosity and 
measurement bias were the parameters adjusted. Previous analyses of the experimental 
data have used graphical techniques to determine vortex-flow parameters. The modified 
Newton-Raphson method of analysis eliminates considerable manual labor and yields more 
consistent vortex-flow parameters. 

The assumed form of vortex flow fit well the measured velocities for the numerous 
sets of data, both flaps up and down for the C-5A airplane. The resulting values of total 
circulation, however, were about two-thirds that expected of a wir^ with an elliptical, load- 
ing. A partial explanation of the less than expected circulation is a dip in the spanwise 
lift distribution at the airplane’s center line. The distance between the trailing vortices 
at the smallest times encountered is somewhat less than that expected for an elliptical 
wing loading. As the time and the distance behind the C-5A airplane increases, the posi- 
tions of the vortex centers become more irregular. Although there is considerable scat- 
ter, the radius of the vortex cores appears to grow approximately as the square root of 
time. The initial core radius, however, is considerably smaller than expected when com- 
pared with values predicted using methods of Spreiter and Sacks and of Betz. The value 
of the effective eddy viscosity which would result in the measured core sizes is higher 
than the values obtained by most other investigators but is consistent with past results 
when the Reynolds number is considered. 



INTRODUCTION 


The hazard of flying in the wake of another airplane, particularly the wake of a large 
airplane, has long been recognized. Usually the problem is avoided by maintainii^ large 
distances between airplanes. In formation flight and during landing and taking off at busy 
airports, however, the distances between airplanes are short so that the vortex flow in 
their wakes can cause dangerous upsets. Unfortunately, efficient use of an airport 
requires airplanes to land and take-off at frequent intervals, the result being that the dis- 
turbing influence of the trailing vortices of airplanes has considerable effect on the safety 
and economic operation of these airplanes, particularly in the terminal area. Conse- 
quently, there is considerable interest in establishing a better understanding of the vortex 
flow in the wake of an airplane so that its disturbing influence can be reduced. For this 
reason the inflight measurements of the vortex flow in the wake of a C-5A heavy transport 
were made (ref. 1) and are being analyzed in this paper. 

Considerable work has been performed in developing mathematical models of vortex 
flow. For example, in reference 2 a finite system of "horseshoe vortices" was used to 
approximate a particular wing loading, which in turn enabled the calculation of the 
"potential flow" at any point. By allowing the trailing vortex filaments to move with the 
induced flow, a "roll-up" of the vortex sheet into two dominant vortices is obtained (ref. 3) 
similar to that observed experimentally. There is, however, some question as to the 
validity of this lumped parameter model as far as the roll up is concerned since, as 
Jordan (ref. 4) has recently pointed out, an elliptical wii^ loading represented by a dis- 
tributed system of vortices does not roll up at all without a disturbance. Also, since the 
viscosity of the air is neglected in these potential flow models, these vortices do not decay 
with age and one must resort to other formulations. 

Lamb, in reference 5, gives a solution for a two dimensional, symmetrical, vortex 
having incompressible, laminar viscous flow. From initial conditions of infinitesimal 
core radius and infinite tangential velocity the solution shows that as time increases, the 
core radius increases, maximum tangential velocity decreases, and total circulation 
remains constant. Time behind the vortex generating airplane can, of course, be 
expressed in terms of distance. While the solution of reference 5 contains certain limi- 
tations imposed by the boundary conditions which do not allow interaction between layers 
of air at different distances behind the generating airplane. Lamb's solution remains a 
useful description of the decay of the vortex flow, particularly at large distances. Super- 
position may be assumed to obtain the total flow for both trailing vortices. Also, it is not 
necessary to use the same initial condition as that used by Lamb since a finite core radius 
and tangential velocity can be used by simply shifting the time reference. A rationale is 
needed, however, for the initial core size. Spreiter and Sacks suggest in reference 6 a 
way of determining a minimum core radius by equating the kinetic energy of the trailing 
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vortices to that created by the induced drag of the airplane. It is also possible to use the 
work of Betz in reference 7 to predict the initial core radius,. But both methods predict 
an initial core radius considerably larger than that observed. 

Because the actual vortex flow of an airplane may be turbulent, the assumption of 
laminar flow in Lamb's solution would not hold. It is common practice to use an empiri- 
cal parameter, the effective eddy viscosity. Values of effective eddy viscosity which have 
been observed experimentally (refs. 8 to 12) range from 7 to 2000 times that of the kine- 
matic viscosity for air. Some investigators (refs. 8 and 13) have suggested that the effec- 
tive eddy viscosity is related to Reynolds number, and the experimental results of several 
investigations (refs. 14 to 18) support this view. It is difficult to obtain accurate values 
of the effective eddy viscosity in wind tunnels, especially at Reynolds numbers approach- 
ir^ flight values, because the test-section size restrictions prevent measurements from 
being made at the proper distances downstream of the model. 

It was because of the hazard of trailing vortices and the limited understanding of 
them that a specially instrumented T-33 airplane was used to repeatedly probe the trail- 
ing vortices of a C-5A airplane. A simplified, graphical analysis of these test results 
have been reported in reference 1. In the present paper the same data is subjected to a 
more sophisticated numerical analysis to determine the characteristics of the vortex flow 
more accurately and to compare these results with analytical models of vortex flow. It is 
intended that an understanding of the phenomenon would be enhanced. 

The basis of the analysis of the vortex flow in this paper involves the automatic 
adjustment of the parameters in a flow equation so that the calculated flow match the mea- 
sured flow. Because of the nonlinear nature of the flow equation, a modified Newton- 
Raphson method (ref. 14) was employed to minimize the mean-square difference between 
the measured and theoretical flow velocities. Having applied this iterative procedure, the 
resulting experimental values, location of vortex center, circulation, tangential velocity, 
core radius, and eddy viscosity could be compared with theory. 

SYMBOLS 


a acceleration, m/sec^ 

b bias in acceleration measurements, m/sec2; wing span, m (see fig. 1) 

b’ distance between vortex centers, m 

c unknown parameter vector; wir^ chord (see fig. 49), m 

c mean aerodynamic chord (see fig. 1), m 
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local lift coefficient 


acceleration due to gravity, m/sec^ 

mean-square error 

iteration index 

roll rate, deg/sec 

pitch rate, deg/sec 

position vector, m 

Reynolds number, 

yaw rate, deg/sec; radial distance, m (see fig. 3) 

radius of vortex core, m (see fig. 3) 
o 

wing area, m* 

translation matrix 

time behind C-5A airplane, sec 

flow velocity relative to inertial axes, m/sec 

velocity of vortex-generating C-5A airplane, m/sec 

measured flow velocity relative to inertial axes, m/sec 

initial measured flow velocity, m/sec 

velocity of probing T- 33 airplane, m/sec 

measured flow velocity relative to probing aircraft, m/s 


components of flow velocity, m/sec 



Vy^O’^z,o initial values of Vy and Vz, m/sec 
^^z,o 

, drift in Vy and Vz, 1/sec 

dy dy 

tangential flow velocity, m/sec 

max ' maximum tangential flow velocity at core boundary, m/sec (see fig. 3) 
x,y,z distance components, m 

a angle of attack, deg 

^ angle of sideslip, deg 

r circulation, mVs^c 

Fq theoretical circulation, mVsec 

B pitch angle, deg 

V kinematic viscosity, mVsec 

effective eddy viscosity, m2/sec 
p density, kg/m^ 

(p bank angle, deg 

heading angle, deg 

Superscript T denotes the transpose; subscripts 1 and 2 denote corresponding 
vortex. 

DESCRIPTION OF FLIGHT TEST 
Airplane Description 

The airplane generating the trailing vortices analyzed in this study was the C-5A, a 
very large transport. A three-view drawing of the C-5A is given in figure 1 and the more 
pertinent dimensions are listed. The mass of the airplane rar^ed from 173 700 kilograms 
to 261 500 kilograms for the subject flight tests. The vortices were visible because of 
entrained engine smoke. 
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Flight Measurements 

The vortex flow behind a C-5A airplane was probed using a specially instrumented 
T-33 airplane which measured linear accelerations, angular rates, flow vane deflections, 
Euler angles, and airspeed. The manner in which the probes were made is shown in fig- 
ure 2. Passes of two types were used, those that were perpendicular to the flight path of 
the C-5A airplane and those that crossed at a shallow ar^le relative to the flight path of 
the C-5A airplane. 

The raw data obtained from the instruments on board the T-33 airplane were 
recorded using oscillograph flight recorders. The data were processed to obtain mea- 
surements of the vortex-flow velocities relative to an inertial reference frame using the 
equation 


Vm = t(Vv + W X r) + |r(a - b) + ^ dt - Vm,o 


where Vm is measured flow velocity and 


R 


Vv = 


-Vp cos a cos ^ 
-Vp sin /3 


\ 


\ 


-Vp sin a 


/ 



vector distance from accelerometer location to vane location 


a acceleration vector at accelerator location 


g gravity vector 

Vm o vector chosen to make V^ =0 at T = 0 

b bias vector chosen to make Vm = 0 at end of run 
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“cos d cos l// 

T = cos Q sin 
_-sin e 

The exact alinement and location of the origin of these inertial axes was different 
for each run or pass of the probing airplane. As figure 2 indicates, the initial position 
attitude and direction of flight of the probing airplane establish the inertial axes used for 
each set of data analyzed. For further details see reference 1. 

Although considerable care was given to insure accurate T-33 instrumentation, it 
was necessary to employ a special filter during data processing to partially compensate 
for the second-order dynamics of the flow vanes and airspeed system. The raw data were 
faired manually to suppress boom bending effects in the vane deflections. 

One additional difficulty is a slight but occasionally important time shift that existed 
between the measured quantities due to differences in the response characteristics of the 
instruments. Although an adjustment was made in the manual analysis of reference 1, the 
adjustment was not made in the analysis using the modified Newton-Raphson method. 

ANALYSIS 

Vortex Flow Equations 

The vortex-flow equation used to fit the flight measurements is based on that given 
by Lamb in reference 5 for two-dimensional, laminar, viscous flow. The tangential veloc- 
ity for a single, isolated vortex is given by 

Vn = — (l - exp—) 

^ 27rr\ 4i/t/ 


The total flow velocity is assumed to be the superposition of two dominant trailing vortices 
with centers at and y 2 »Z 2 - The three components of the calculated flow velocity 

are given by 



sin 4> sin 0 cos \{/ - cos <f> sin \{/ cos <f) sin 6 cos i// + sin </> sin 

sin <p sin 6 sin \p + cos (f) cos xl/ cos 4> sin 9 sin ip - sin 0 cos \p 

sin <p cos 6 cos <p cos 6 
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where 


^1 






The terms V , V , - and — ^ are called bias terms and are added to 

y,o z,o dy dy 

account for uncertainties of the inertial reference frame. 

There are, of course, many conditions which pertain to Lamb's solution which are 
not met. If the flow were actually symmetrical and laminar, then the value of the kine- 
matic viscosity u would simply be that for air. Unfortunately, the vortex flow in the 
wake of an airplane may be turbulent, and the initial roll-up process which forms a pair 
of almost symmetrical vortices may not be complete. Consequently, a value of that 
is many times that for air is often required to match the observed radius of the vortex 
core (ref. 11). This empirical value of i' is termed the effective eddy viscosity 
The assumption of two-dimensional flow ignores the effect of shear between layers of air 
at different distances from the vortex-generating airplane and does not account for the 
axial flow Vx- Superposition is assumed to be valid for the flow of the two vortices 
although it is known that the Navier-Stokes equations are nonlinear. 

McCormick in reference 16 has advanced an alternate form of flow equation which 
has been motivated by the empirical models of turbulent boundary layers and is compared 
with Lamb’s solution and Spreiter-Sacks approximation in figure 3. It is clear that 
McCormick’s formulation does not behave properly as the radius approaches zero. The 
limited information available, however, precludes a clear decision as to which flow equa- 
tion might be most accurate. Lamb’s solution was chosen simply because of its wide use 
and its analytical rather than empirical basis. 


Modified Newton-Raphson Method 

The purpose of the analysis is to determine the parameters r and in the flow 
equation so that the calculated flow is as close as possible to that measured. It is nec- 
essary also to determine the location of the centers of the two vortices and four bias 
terms to account for uncertainties in the reference frame, making a total of 10 unknown 
parameters: 


T Q Q 
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It is apparent from the nonlinear form of the flow equation that an iterative pro- 
cedure is required. A modified Newton-Raphson method (ref, 19) was selected because 
of its successful application in system identification problems. The criterion used in fit- 
ting the calculated flow to that measured is the mean-square error: 

N 

J [Vm(‘n) - V(tn)]'^[Vni(t„) - V(t„)] 
n=l 

Estimates of the unknown parameters which form the vector c can be found by 
Solving for the vector c which minimizes J. Since this condition is met when the 
gradient of J with respect to the unknown parameter vector c is zero, then c will 
be a root of the equation 

VcJ-0 

If the gradient is taken to be equal to 

= ^c^k + ^c^Jk(^k+l ■ *^k) ^ 
then solving for gives 

r 2 i‘‘ 

*^k+l " *^k ■ -^kj ^cJk 

This is an application of the Newton-Raphson method. 

For the problem considered, the first gradient is 

N 

V^J = -2 ^ V(.v'^(vm - v) 

n=l 

and the second gradient is 

N N 

= -2 ^ - 2 ^ Ve2v(v^ - v) 

n=l n=l 
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where the second term is negligible. In the interest of computational efficiency and pro- 
graming simplicity, the second term in the expression for the second gradient is neglected. 
The recursive formula that is consequently used is 


‘^k+l " “ 



DISCUSSION OF RESULTS 


Comparison of Methods of Analysis 

An example result of applying the Newton-Raphson method to obtain estimates of the 

unknown parameters y^ 73 , ^ 2 ’ 17 ’ T^’ ^y* ’ ^z,o’ ^y’ 

shown in figure 4. Shown are the mean-square error (dashed) and the 10 unknown param- 
eters as a function of the iteration index. It is interesting that for the example shown, the 
mean- square error would appear to converge after seven iterations but the bias term 
V,, n required 20 iterations to settle completely. In general, the closer the probing air- 
plane passes to the center of the vortices, the greater the reliability of the estimates and 
the fewer the number of iterations required. It is reasonable to expect this because of 
the insensitivity of the mean-square error to changes in the unknown parameters that 
occur at the large distances from the vortex centers. Consequently only runs which 
reached the edge of the core were used to determine the parameter 

An example comparison of the measured and theoretical vortex flow that results 
from applying the Newton-Raphson method is shown in figure 5. The data are for a per- 
pendicular probe of the wake of a C-5A airplane with flaps up and a vortex age of about 
24 sec. The theoretical values match those measured fairly well, considering both com- 
ponents must be matched simultaneously. By noting the sign of the lateral velocity at the 
points at which the vertical velocity is zero, it can be seen that the probing airplane passed 
above the center of the first vortex and below the center of the second vortex. Because the 
probing airplane happened to pass almost as close to the center of the second vortex as the 
first, the velocity functions are almost symmetrical. 


At the bottom of figure 5 is a plot of nondimensional tangential velocity against non- 
dimensional distance. The distance from the core has been divided by ^4i^ and tangen- 
tial velocity by rQ/j^27r^4i^^, The curve without symbols represents the calculated veloc- 


ity; the curve with symbols, the measured. Both curves start on the right side, cross the 
origin as the probing airplane passes over or under the center of the first vortex center. 
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then double back over the origin as the second vortex is passed. The plot is useful in 
assessing whether a core penetration has been achieved by noting the minimuna distance 
achieved. 

Convergence is not always achieved in applying the Newton-Raphson method to a 
problem that is very nonlinear in the unknown parameters. If the starting values of the 
parameter differ greatly from those which minimize the mean-square error J, the itera- 
tive process will not converge. Preliminary work that is reported in reference 20 indi- 
cates the starting values of the vortex center positions need to be within 7 m of the final 
values. Convergence has not been difficult, however, to achieve in practice. One inter- 
esting feature of this particular problem is that the gradient of J approaches zero not 
only as convergence is obtained but also if the iterative process diverges. This is because 
of the insensitivity of the mean-square error to parameter changes at large distances from 
the centers of the vortices. One is assured of a legitimate solution however, if the cost 
subsides, if parameter values stabilize, and if the gradient of this cost continues to 
approach zero. 

A graphical method, which was exnployed in reference 1 and is illustrated by fig- 
ure 6, can be used to determine vortex centers. Lines are drawn perpendicular to the 
resultant flow. In the case of an isolated symmetrical vortex, these lines should pass 
through the center of the vortex. Errors in the flow measurements and the presence of 
at least two vortices limit, however, the accuracy with which the vortex centers can be 
determined. The motivation for applying the modified Newton-Raphson method to analyz- 
ing vortex-flow measurements was to eliminate the manual operations involved in the 
graphical approach and to improve the accuracy of the resulting estimates of core radii 
and effective eddy viscosity. 

The flow velocities which result from applying the two methods of analysis are com- 
pared with that measured in figures 7(a) and 7(b). In the case of the graphical method, 
the values of core radius and tangential velocity given in reference 1 were used with 
Lamb’s equation to compute the velocities shown. The calculated plots of the flow veloc- 
ities for the modified Newton-Raphson technique are a best fit of the flight measurements. 
The positions uf the vortex centers, eddy viscosity, circulation, and bias terms on the 
measured flow velocities were adjusted by the algorithm to minimize the mean-square 
error. As expected, the fit of the measured flow velocities is better In the case of the 
modified Newton-Raphson method, simply because that is the purpose of its use. In the 
case of figure 7(a) the difference is small, whereas in figure 7(b) the fit produced by the 
Newton-Raphson method is significantly better. Disagreement between the calculated and 
measured functions can be due to the limitations of using Lamb's solution as well as mea- 
surement errors. 
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Comparison of Calculated and Measured Flow 

Figures 8 to 44 show the calculated vortex flow usii^ a modified Newton- Raphson 
method compared with measured flow. Figures 8 to 22 are for the flaps up on the C-5A 
airplane, figures 23 to 44 are for flaps down. The figures axe in order of increasing time 
behind the C-5A airplane. Table I contains a list of the vortex and generating airplane 
parameters that correspond to the figures. 

The comparison between the calculated and measured flow is good in most cases, 
especially considering that both the lateral and vertical velocity components are fit simul- 
taneously. Since the data for each run were taken for less than 1 sec, adequate instru- 
mentation was a concern. A periodiclike error seen in many of the cases is attributed to 
nose boom bending. This error can be noted in figure 9, for example. 

Some of the cases for flaps down indicate additional peaks in the flow velocity so 
that instead of two dominant trailing vortices there are four. This is most evident in 
figure 44. Unfortunately, the boom bending causes peaks at about the same period so 
that only speculation can be made. In the Newton-Raphson analysis, however, no provi- 
sion for boom bending was made and only two vortices were assumed; as a result, the 
fit of the measured data was poor. 

VORTEX CORE RADIUS AND EFFECTIVE EDDY VISCOSITY 

It is usually assumed that the rollup of a vortex sheet behind a wing is essentially 
complete after several seconds. Further changes in the radius of the core of the pair of 
trailing vortices have then been attributed to viscosity of the air. For two-dimensional, 
symmetrical, viscous flow. Lamb shows the core radius r^ to be directly related to the 
kinematic viscosity by the expression 

r^ = 2.24 

Unfortunately, the core radii calculated using this expression have not agreed with obser- 
vations of vortex flow. Because of this, the practice has been to replace the kinematic 
viscosity with an empirical parameter, the effective eddy viscosity, and to adjust its 
value until the calculated core radius agrees with that observed. If this practice is 
applied to the results of the analysis presented herein, the comparisons shown in fig- 
ures 45 to 50 can be obtained. 

Since core radius, circulati.on, and maximum velocity are interrelated, it was nec-. 
essary to consider the fairings of figures 45, 48, and 50 simultaneously to ensure the con- 
sistency of the faired results. In fairii^ the core-radius results of figure 45, it was 
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assumed that the rate of growth given by Lamb’s solution would hold but that a finite, 
rather than an infinitesimal, core radius would exist at the time zero. It seemed reason- 
able to use the same value of effective eddy viscosity but not the same initial core radius 
for both the flaps-up and flaps-down cases. The equation of the fairings are given in fig- 
ures 45(a) and 45(b). 

Attempts were made to calculate the initial core radii using both the methods of 
Spreiter and Sacks of reference 6 and the radius of gyration of the vortex system dis- 
cussed by Betz in reference 7. Both methods indicated the radius to be about 7 m to 8 m, 
about three times as large as that observed at points closest to the C-5A airplane. 

Because of this discrepancy an empirical value was used which best fit the observations. 

The overall value of effective eddy viscosity which corresponds to the fairir^s shown 
on figures 45(a) and 45(b) is 634 times the kinematic viscosity u of air. Figure 46 shows 
that this ratio of is considerably greater than the corresponding experimental val- 

ues taken from references 8, 9, 10, and 12 but is comparable to that observed by Rose and 
Dee in reference 11. Squire in reference 13 has suggested that the ratio of v^jv should 
be proportional to a form of Reynolds number, R^ = while Owen argues v^jv 

should vary proportional to the square root of Rjj. Figure 46 shows the results of sev- 
eral investigations on a plot of v^jv versus Rj^ and indicates a definite trend with 
Reynolds number. It would appear that being proportional to the square root of 

as suggested by Owen, is the better description of the trend shown. A line with slope 
one-half has been drawn through the value of v^jv obtained in the subject work. As 
is decreased, it would be expected that the slope will become zero and will approach 

one. 

It is questionable that the entire growth of the vortex core can be attributed to vis- 
cous effects, even though it is common practice. It is disturbing that a ratio of v^jv on 
the order of 600 is needed to match the observed core size; it is also disturbing that the 
core does not quite obtain the size expected from potential flow considerations. It is 
suspected that at least part of the observed growth in the core radius is due to the original 
rollup of the vortex sheet and is not entirely due to viscosity. Unfortunately, the informa- 
tion available is insufficient to resolve this question. 

Distance Between Vortex Centers 

The distance separating the trailing vortex centers as determined from the Newton- 
Raphson analysis is shown in figure 47. The value calculated using an elliptical spanwise 
lift distribution, which is 77/4 times the span, is somewhat greater than that observed. 

As the distance and time behind the C-5A airplane is increased, a dispersion occurs which 
caused considerable variance in the distances measured. This is thought to be due to dis- 
turbances in the atmosphere and due to inherent instabilities of the pair of trailing vortices. 
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Parks (ref. 21) and many others have modeled and analyzed such phenomena and have 
shown that the trailing vortices snake in planes at 45° to the vertical, with a wave length 
of 7.18 times the distance initially separating the vortices. Parks also gives a rate of 
divergence of the snaking. Asymptotes have been drawn on figure 47 to show the rate at 
which the snaking motion would grow from the same initial amplitude. It is interesting 
that the divergence envelopes resemble the growth of the dispersion of the distance 
between vortex centers. In smooth air the snaking results in the formation of "doughnuts" 
or rii^ vortices. If the wavelength of the snaking is taken as the half of the circumference 
of a circular ring, the maximum distance could reach 3.6 times the span. It can be seen 
in figure 47 that this distance is considerably greater than the largest distance measured. 

Total Circulation 

Another of the results of the analysis was a measure of the total circulation for each 
of the trailing vortices. The values obtained have been normalized using the value 
expected for an elliptical wing loading and are plotted in figure 48 for both flaps up and 
flaps down. Note that the total circulation obtained is generally less than that expected 
for an elliptical wing loading. Similar results have been obtained by many other investi- 
gators. It is believed that one reason lies with the fact that the spanwise lift distributions 
involved are not elliptical but have a dip at the plane of symmetry. Since the total circu- 
lation depends only on the center-line local lift coefficient by the relationship 

c,c 

To = j- pVg 

a dip in the wing loading distribution at the center line would produce less circulation than 
for an elliptical distribution. In the case of the C-5A flight results, a dip is known to 
occur for the wing and is expected to be amplified by the fuselage and tail. Figure 49 
shows the wii^ spanwise lift distributions furnished by the Lockheed-Georgia Company. 
The value of the local Uft coefficient at the plane of symmetry that would be necessary to 
agree with measured circulation is indicated by a circle for comparison. Because of 
interference of the stii^ mount used in the tunnel tests by Lockheed-Georgia, it was not 
possible to get complete data at the plane of symmetry. 

The peak flow velocities divided by the circulation of an elliptical wing loadii^ is 
presented for the C-5A airplane in figure 50. Vertical leaders on the points indicate the 
measured maximum values while the solid symbols give the maximum calculated velocity 
for each run. The open symbols are also the maximum calculated velocities; but since 
core penetration was not achieved, there is less confidence in their values. Except for a 
few points having considerable discrepancy the data is quite consistent. The fairings used 
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observed the direct relationship between the circulation of figure 48 and the core radius 
of figure 45. 


CONCLUDING REMARKS 

A modified Newton-Raphson method has been used successfully to determine esti- 
mates of the total circulation, effective eddy viscosity, and the locations of the centers of 
the trailing vortices of the C-5A airplane based on flight measurements. Lamb’s solu- 
tion has been used to model the vortex flow, and the comparison between the calculated 
flow and measured flow is good. The convergence properties when applying the modified 
Newton-Raphson method to such a nonlinear problem were good but required up to 20 
iterations. 

The vortex-core-radius values that were obtained using the analysis had consider- 
able scatter but appear to grow in a manner similar to Lamb’s solution but with a finite 
core radius initially. The initial core size obtained using the method of Spreiter and 
Sacks is considerably larger than that measured. If the growth in core size is attributed 
solely to viscous effects, the value of effective eddy viscosity obtained is greater than the 
results of most other investigations but follows a trend indicating effective eddy viscosity 
to be proportional to the square root of Reynolds number. The distance separating the 
vortices at the smallest times encountered is somewhat less than that predicted for an 
elliptical wing loadii^. At greater distances behind the C-5A airplane, dispersion exists. 
It is suggested that the dispersion is due in part to the unstable snaking action modeled by 
Parks and others. The total circulation obtained is less than that expected for an ellipti- 
cal spanwise lift distribution and decreases with time and distance behind the C-5A 
airplane. 

Previous analyses of the experimental data have used graphical techniques to deter- 
mine vortex flow parameters. The modified Newton-Raphson analysis reduces consider- 
ably the manual effort of graphical techniques and is believed to give more accurate 
results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 3, 1974. 
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APPENDIX 


GRADIENT EXPRESSIONS 

The formulation used for the gradient of the flow velocity V with respect to the 
unknown parameter vector c is as follows: 
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APPENDIX - Concluded 
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TABLE I.- VORTEX AND GENERATING AIRPLANE PARAMETERS 


Symbols used in 
figs. 45, 47, 
48, and 50 
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97 
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^Filled symbols denote core penetration. 

Adjusted to be consistent with distance between vortex centers, 
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ng of the C-5A heavy jet transport. 
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Figure 4.- Convergence properties of modified Newton-Raphson method. 


,00 


(a) Velocity versus distance. 
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(b) Nondimensional velocity versus nondimensional distance. 
Figure 5.- Comparison of measured and theoretical vortex flow. 
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Figure 7.- Comparison of methods used to analyze vortex-flow measurements 

for two examples. 
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Figure 8.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex age of 12 seconds. 
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Figure 9. - Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex age of 19 seconds. 
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Figure 10.- Fit of measured data using modified Newton-Raphson method 
for flaps up and vortex age of 24 seconds. 
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Figure 12.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up (slide through) and vortex age of 32 seconds. 
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Figure 13.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex age of 34 seconds. 






Figure 14.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up (slide through) and vortex age of 42 seconds. 






Figure 15.- Comparison of Newton-Raphson fit with measured flow velocities 


for flaps up (slide through) and vortex age of 46 seconds. 
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Figure 16.- Comparison of Newton-Raphson fit with measured flow velocities 


for flaps up and vortex age of 60 seconds. 
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Figure 19.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex a^e of 70 seconds. 
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Figure 21.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex age of 104 seconds. 
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Figure 22.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps up and vortex age of 114 seconds. 
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Figure 23.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 22 seconds. 
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Figure 24.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 30 seconds. 
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Figure 25.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 32 seconds. 
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Figure 27.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 41 seconds. See table I for variable 
changes. 
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Figure 28.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 41 seconds. See table I for variable 
changes. 
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Figure 29.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 42 seconds. 




Figure 30.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 44 seconds. 
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Figure 31.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 59 seconds. 
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Figure 32.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 60 seconds. See table I for variable 
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measured flow velocities 
} table I for variable 
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Figure 35.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex £^e of 64 seconds. 
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Figure 36.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex of 68 seconds. 






Figure 38.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 82 seconds. 
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Figure 39.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 85 seconds. 
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Figure 40.- Comparison of Newton-Raphson fit with measured flow velocities 


for flaps down and vortex age of 92 seconds. 
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Figure 41.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 109 seconds. 
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Figure 42.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 115 seconds. 
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Figure 43.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 116 seconds. 






Figure 44.- Comparison of Newton-Raphson fit with measured flow velocities 
for flaps down and vortex age of 133 seconds. 
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Figure 45.- Core radius as a function of time. (See table I for symbol notation.) 
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Figure 47.- Distance separating vortex centers. (See table I for symbol notation.) 
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Figure 48.- Normalized circulation versus time. (See table I for symbol notation.) 
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Figure 50.- Maximum-circumferential-velocity circulation ratio versus time. 

(See table I for symbol notation.) 
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